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ABSTRACT 

The possibility that the stellar initial mass function (IMF) arises mostly from cloud 
structure is investigated with fractal Brownian motion (fBm) clouds that have power- 
law power spectra. An fBm cloud with a realistic projected power spectrum slope of 
(3 = 2.8 is found to have a mass function for clumps exceeding a threshold density 
that is a power-law with a slope of a = 2.35, the same as in the Salpeter IMF. Any 
hierarchically structured cloud has a clump mass function with about the same slope. 
This result implies that turbulent interstellar clouds produce dense substructure with 
the observed pre-stellar core mass function built in from the start. Details of the clump 
formation processes are not critical. The conversion of clumps into stars involves a 
second step. A one-to-one correspondence between clump mass and star mass is not 
necessary to convert the clump mass spectrum into an IMF with the same power- 
law slope. As long as clumps have an internal stellar IMF from sub-fragmentation, 
protostellar accretion, coalescence and other processes, and the characteristic mass for 
this internal IMF scales with the clump mass, then the IMF slope above the minimum 
characteristic mass will equal the clump mass slope. A detailed review of IMF models 
illustrates the prominence of cloud structure as a major component in a wide class of 
theories. Tests are proposed to determine the relative importance of cloud structure 
and competitive accretion in the IMF. 
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for giant molecular clouds, molecular cloud clumps, star 
clusters, and individual stars, we would like to understand 
the relationship between the structure that is viewed as 
a featureless power spectrum and the structure that is 
inferred to be a collection of discrete objects. 



IStiitzki et al.l l| 19981 ) suggested that for Gaussian- 
shaped clouds with a mass-radius relation M oc B? and 
a power law mass function of slope —a, the slope — /3 of 
the power spectrum is given by f3 = 7 (3 — a) . This re- 
lation comes from an integral over the cloud mass func- 
tion, dn oc M~ a dM, of the contribution to the power spec- 
trum, dP(k), from a cloud of mass M, which is dP(k) oc 
M 2 exp(— 0.5k 2 R 2 )dn. The result is the frequency (k) de- 
pendent part of the integral over M 2 ~ a exp(-0.5k 2 R 2 )dM 
which is fc - ' 3 for the above value of fi. IStiitzki et al.l (|l998l ) 
noted that the value of /3 ~ 2.8 they observed corresponded 
to typical a ~ 1.6 and 7 ~ 2 for molecular clouds. Observa- 
tions of /3 for various regions of the Milky Way are in Table 
1; generally /3 ~ 2.8. 



1 INTRODUCTION 

1.1 Clump Mass Functions and Power Spectra 

Interstellar gas emission has a power-law power 
spectrum su ggesting structure on a wide range 



of scales (ICrovisier fc Dickevl Il983l: fCreer] [l993: 



Dickey. McCiure -Grifnths, Stanimirovic, Gaensl er fc Greenl 



200l| ; IMiville-Deschenes, Joncas. Falgarone fc Boulangerl 
2003: IMiville-Deschenes. Levrier fc Falgarone! 12003" : 



Miville-Deschenes. Martin fc et ail |2010T ). Cloud-fitting 



algorithms that are applied to this structure find 
power law mass functions for the emission peaks, which 
are usually identified as clouds, clumps, and core s 
l|Williams. de Geus fc Bht3ll994l ; IStiitzki fc Giiestenlll99ch . 
Because there is considerable interest in the origin of 
various interstellar and stellar mass functions, such as those 
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Table 1. Observations of Power Spectrum Slope ft for various regions of the Milky Way. 



Region 


Type of Observation 


Power Spectrum Slope 

(ft) 


Reference 


Foreground of C as A 


HI 21 cm absorption 


2.75 ± 0.25 


Deshpande et al. f 20001 






2.86 ±0.1 


Roy et al. (2010) 


Perseus, Taurus, Rosetta clouds 


12 CO 


2.74 ±0.08 


Padoan et al. (20041 


Perseus cloud 


13 co 


2.86 ±0.1 


Padoan et al. f 20061 




12 CO and 13 CO 


w 3.1 


Sun et al. C20061 


Perseus spiral arm 


HI 21 cm 


2.2 to 3.0 


Green (1993) 


Ursa Major high-latitude cirrus 


HI 21 cm 


3.6 ± 0.2 


Miville-Deschenes et al. (2003) 


Polaris Flare 


12 CO 


~ 2.8 


Stutzki et al. ( 19981 




FIR 


2.7 ±0.1 


Miville-Deschenes et al. (2010) 


Several molecular clouds 


12 CO and 13 CO 


2.5 to 2.8 


Bensch et al. (20011 




100 /im 


2.9 to 3.2 


Gautier et al. ( 19921 


The Fourth Galactic Quadrant 


HI 21 cm 


^ 4 


Dickev et al. (20011 


The Gum nebula 


8, 24, and 70 ^m 


2.6 to 3.5 


Inealls et al. (20041 



iHennebelle fc Chabrierl (|2008l ) derived the cloud mass 
distribution from a power spectrum in a different way. They 
assumed that the logarithm of the density, rather than 
the density, has a po wer-law power spectrum (following 
iBeresnvak et al] (|2005l l). and they considered a log-normal 
density pdf, V(p,R), for scale R. Then they equated the 
total mass of clouds more massive than M from the inte- 
gral over the mass function, J M'n(M')dM' , to the to- 
tal mass in the pdf at densities higher than p — M/R 3 , 
which is M pdi (p,R) = J V(p' , R)dp' . This gives n(M) = 

M- 1 dM pdt (p,R)/dM = M- 1 (dR/dM)dM pdf (p,R)/dR. 
The pdf V depends on R because the dispersion a of 
the log normal function was assumed to depend on R, as 
a 2 = a 2 (l - (R/L) n '- 3 ^ for n' > 3. That is, the den- 
sity fluctuations become smaller when averaged over larger 
scales, up to the system scale L. This dependence comes 
from the assumption a 2 oc f^!^ k~ n d 3 k where n' is the 
(negative) slope of the power-law power spectrum for log(p). 
The result is that n(M) oc M~ 1 a~ 

M -2 M n'/3-l = M - a for R << L 

This is about the same as in IStutzki et ail Jl99ct). who had 
q = 3 — /3/3 for M « fi 3 , considering that ft in Stutzki et al.l 
was the power spe c trum slope for density while n' 



3 {dR/dM){da 2 /dR) oc 
, where a = 3 — n'/3. 



Hennebelle k, Chabrierl ((20 08) was the power spectrum 



slope for log-density. 



The IHennebelle fc Chabrierl |200Sl ) model has a differ- 
ent relation between a and ft when n' < 3, because then 

a 2 — a 2 ^(L/R) 3 ~ n — lj diverges as R goes to zero. When 

n' > 3 as in the previous paragraph, a approached a con- 
stant for small R and its presence as a~ 3 in the expres- 
sion for n(M) contributed no additional mass dependence; 
only da 2 j dR contributed, making n(M) oc M~ 2 Rda 2 j dR oc 
M~ 3+n / 3 as above. When n' < 3, however, a contains a 
strong M dependence and the a~ 3 term is important, as is 
the tj 2 in another term, (8 C + a 2 /2), which comes from the 
log-normal density distribution for threshold density S c . This 
makes n(M) oc M^a' 1 Rda 2 /dR oc M~ 2 a oc M- 2 ' 5+n ' /6 . 



1.2 Clump Mass Functions and Hierarchical 
Structure 

For hierarchically structured objects, the mass function 
n(M) always has a power law slope a = 2 if the total mass 
of clouds in a logarithmic mass interval around a partic- 
ular mass is the same for all masses, i.e., for all levels in 
the hierarchy. Then MdN(M)/dlogM — constant for mass 
function N in logM intervals, giving N(M) oc M~ , or 
n(M) oc M~ 2 for linear M intervals. This a — 2 result 
can be derived in many ways. Clouds with an average of 
TV substructures per structure in all levels L of the hierar- 
chy have a number of objects per logarithmic mass inter- 
val that increases as N L and a mass per ob ject that de - 
creascs as 1/N L , making dN/dlogM oc M _1 |Flecklll996l ). 
Similarly, the probability of selecting out of all log M lev- 
els a particular structure with mass M is proportional to 
the number iV of these objects, which is also oc M for 
logM intervals. The same result applies to packing scale- 
free structures in a volume. For wavenumber k proportional 
to the inverse of the size, the number of structures be- 
tween k and k + dk that fit into a space with dimension 
D is £(k)dk oc k D ~ 1 dk, and the mass of each is M oc k~ D . 
Converting wavenumber counts into mass counts using the 
one-to-one correspondence, ^(k)dk — n(M)dM, we derive 
n(M) oc k u ~ x dk ldM oc k 2D oc M~ 2 , indepen dent of di- 
mension D. The IHennebelle fc Chabrieil (|2008) result for 
non-gravitating clouds is shallower than M~ 2 by the index 
ra'/3 — 1 = 2/9 for n' = 11/3 because the total cloud mass 
in each log M interval is not constant with M, but decreases 
with increasing M as the pdf becomes narrower and the in- 
tegral of the pdf above a threshold density becomes smaller. 



1.3 The Stellar Initial Mass Function 

The initial mass function for stars (IMF) seems consider- 
ably more difficult to understand than the mass function for 
non-gravitating clouds because the IMF involves the whole 
star formation process, including a wide range of densities 
and diverse physical effects. Still, the IMF in clusters is of- 
ten observed to have a power law component spanning 1.5 
to 2 orders of magnitude in the range from ~ 0.5 Mq to 
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50 Mq. How much of this power law comes from a power- 
law in cloud structure and how much comes from other pro- 
cesses is not currently known. Current explanations for the 
IMF consider cloud fragmentation driven by supersonic tur- 
bulence and self-gravity, clump coalescence, and protostellar 
accretion (see reviews in iBonnell et alll2007l : iDib et al.ll20ld : 
iBastian et al]|2010h . 

IPadoan fc Nordlundl J2002T ) considered an IMF model 
based on turbulent fragmentation of a cloud. They proposed 
that turbulence compresses the gas into layers of thickness 
L oc V _1 for ambient turbulent speed Vq oc LJj, initial 
length Lo, and Kolmogorov exponent jj ~ 1/3. They also 
assumed that the number of stars scales with LjT 3 d log Lo 
from close packing (as above, where we wrote £(k)dk oc 
fc 3 dlogfc for k = 1/L in 3D). Then with M oc pL 3 = 
PqLqL 2 for cubic regions in the compressed layer and den- 
sity p = poLo/L for uncompressed density po, it follows that 
M oc Lq (L/L ) 2 oc Lo" 2 ". Thus dN(M)/dlogM oc L^ 3 oc 
M 3/(2r,-S) ^ M' 1 - 3 when 77 = 1/3, which is close to the 
Salpeter IMF, dN/dlogM oc M" 1 ' 35 . The final IMF was 
assumed to be this clump mass function multiplied by the 
probability that a clump of mass M exceeds the thermal 
Jeans mass; the pdf for thermal Jeans masses comes from 
the pdf for density at a constant temperature. Their clump 
model steepens the mass function slope above the value of 
— 1 for mass-conserving fragmentation because they limit 
the number of stars that can form in each layer of dimen- 
sion L x Lo x Lo to a constant, even though the available 
space scales with Lq, i.e., (Lo/L) 2 cubical volumes fit in this 
layer. If each of these cubes were to form a star, then N(M) 
would be larger by the factor (Lo/L) 2 and we would retrieve 
the purely hierarchical result, dN(M) /dlog M oc M~ x . 

IPadoan et all (|l997T l and lHennebelle fc Chabrier] (|2008T ) 
proposed a different model for cloud fragmentation in 
the self-gravitating case. They start with the expression 
written above for a non self-gravitating cloud, n(M) oc 
M~ 1 dM pc \t (p, R)/dM, but then consider the unstable re- 
gions with average density p larger than the density at which 
the Jeans length equals R. Instead of summing the cloud 
masses above p-R 3 , they sum the cloud masses smaller than 
pi? 3 , because the cloud's mass will be less than the Jeans 
mass at length R when the density exceeds the Jeans-length 
density. The result is n(M) oc M " 1 (dp/dM)P v< j t (p, R) 
where p is the density giving M = M.i. IPadoan et al l [|l997f l 
assume the masses are thermal Jean s mass es, in which case 
M oc p~ 5 . iHennebelle fc Chabrieil |2008T l consider a tur- 
bulent medium, for which Mj ~ GV 2 R oc R 2r,+1 when 
V oc R v , and then p ~ Mj/R 3 oc R 2r, ~ 2 oc M (^-2)/(2 V +i) _ 
In both cases, the mass function is the product of a power 
law and a log-normal. In part of the mass range, the slope 
has the Salpeter value. 

IDib et all |2010r i presented a semi-analytical model 
to describe the simultaneous evolution of the dense core 
mass function and the IMF, considering accretion and feed- 
back by winds. The basic core mass function came from 
IPadoan fc Nordlundl l|2002h . The results showed mass func- 
tions that varied throughout the cloud and changed with 
time. The final functions were those at the time when 
pre-stellar winds disrupted the cloud, and they compared 
well with the core and stellar mass functions in Orion. 



Because the underlying m odel is the core function in 
IPadoan fc Nordlundl (2002), geometric effects play an im- 
portant role in their IMF. 

All of the analytical models just described have an 
q = 2 power law as an underlying mass function, modi- 
fied in the various cases by more or less total mass being 
included in logarithmic intervals as a function of the cloud 
mass. The models that begin with an equation Mn(M) = 
dM p di(p, R)/dM for M p df from the density pdf effectively 
assume that an increment toward larger or smaller mass 
in the mass function corresponds to an increment in den- 
sity. For the IMF models, this means that massive stars 
correspond to low density clouds or parts of clouds. This 
is contrary to the mass segregation that is often observed 
in young clusters, and requires accretion or coagulation as 
a second step. The same mass-density relation results from 
any isothermal model that identifies star mass with Jeans 
mass along the power law part of the IMF. To avoid this 
contradiction, the pr e-stellar gas temperatu re has to increase 
with core mass (e.g. , iKrumholz et afj|201Ch . 



1.4 The Press and Schechter Mass Function 

The cosmologi cal equivalent o f these derivations was con- 
sidered by IPress fc Schechterl (|l974T ). They assumed in- 
finitesimal perturbations in an expanding universe and de- 
termined the mass function of bound objects as a func- 
tion of time. The initial density was assumed to have a 
power spectrum, P(k), and the distribution function, V, 
for density perturbation amplitude 5 inside a given scale 
R was assumed to be Gaussian in 8 with dispersion a 
given by a 2 = (2tt)~ 3 / d 3 kP(k) W 2 (R, fc). The window 
function W(R, k) was introduced by iPadmanabhar] (|l993h 
to avoid a divergence at high k when f3 < 3; a Gaussian 
form W ~ exp (~0.5k 2 R 2 ) was assumed. If 5 C is the crit- 
ical density at the initial time that makes a perturbation 
bound at the present time, then the fraction of bound ob- 
jects at the present time with masses larger than M is 
F(M) = TdS. The co-moving mass distribution func- 
tion is then n(M) = (p /M)dF/dM . This procedure is sim- 
ilar to that followed by IHennebelle fc Chabrier] (|2008h . ex- 
cept that Hennebelle & Chabrier (2008) also have a lower 
limit to fc in the integral for a. This lower limit intro- 
duces a constant part in the expression for a 2 , and causes 
the mass-dependence of a to be negligible when n' (j3 
in our notation) exceeds 3 in the Hennebelle & Chabrier 
(2008) model. Without a lower limit on fc, the cosmological 
mass function is n(M) oc M~ a~ exp (-<5 2 /2ct 2 ), which 
is the Press & Schechter (1974) result. Note the purely 
geometric component, M~ 2 , as in the other mass func- 
tion models discussed above. For initial power spectrum 
P(k) oc k~P and mass M oc fc -3 , the dispersion in the den- 
sity distribution function is a oc M (/3_3)/6 . Then n(M) oc 
M- 2 M (3 - /3)/6 exp(-0.5[M/M tr ] (3 "' 3)/3 ). This mass func- 
tion is a power law truncated by an exponential. The trun- 
cation mass, M tr , is the mass that is at the critical den- 
sity for boundedness at the present time. The slope of 
the power law part, a — 1.5 + /3/6, has a positi v e cor - 
relation between a and /3, unlike the IStiitzki et al.l (|l998l ) 
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and lHennebelle fc Chabrierl l|2008l ) slopes. We return to this 
point again in Section [3721 



1.5 The IMF in Numerical Simulations: 

Competitive Accretion or Cloud Structure? 

Numerical simulations of star formation have a different ap- 
proach. They explain the IMF by analogy to the mass dis- 
tribution of sink particles formed in the s i mulation (e.g. , 



Tillev fc Pudrit j 12007); iPadoan et alj 12007"; iLi et al 



Nakamura . 



Lj |2005l . 



2007; Martcl ct al. 



2006 



12004 
iBatdboogh . 



There are many physical processes involved, and it is not 
generally clear which dominate the IMF. Moreover, the sim- 
ulations differ from each other in the presence of magnetic 
forces, magnetic diffusion, turbulence driving, initial cloud 
boundedness, gas cooling, stellar outflows, starlight heat- 
ing, and other properties, even though they get about the 
same IMF. One thing they have in common is a turbulent 
gas, which means hierarchical structure with an M~ 2 un- 
derpinning for mass functions. But there are other things in 
common too. 

It is instructive to see what differs in the few simulations 
that gave an other IMF. Tw o produced relatively flat IMFs: 
one model in iKlessenl j200lf) that had small-scale turbulence 
driving, and models in lClark et al.l (|2008T ) that had transient 
unbound clouds. What they had in common was isolated 
star formation in non-hierarchical clouds, as determined by 
a short driving scale in the first case and by rapid cloud 
expansion in the second. They also differed from the usual 
simulations in having no competition among protostars for 
the gas: one sink particle formed in each fragment with an 
overall efficiency that was low. 

Competi tive accretion JZinnecker] 

19821; iBonnell. Bate. Clarke fc Pringld l200lt 

Bonnell. Clarke. Bate fc Pringld 1200 ll ) involves cloud 



mass (and ultimately sink particle mass) that comes from 
all over a cloud, indepe ndent of the initia l mass of the clump 
that is accreting (e.g.. Iciark fc Bonnelll [20061 ; ISmith et al.l 
2008). The final star mass then depends on the accretion 
rate integrated over time. This mass can be small if the star 
is ejected from its ga s reservoir early l|Reipurth fc Clarkd 
l200ll ; iBate et al.ll2002| ). or if the accretion is slowed by any 
of a number of pro cesses, such as tidal forces in a cluster 
(Bon neli et alj |2008). low local densiti es, or relative mot ions 
between the protostar and the gas (|Clark et alj 1200*91 ). If 
competitive accretion is dominant over geometric effects, 
then it should give the Salpeter IMF in a uniform static 
cloud that is artificially seeded with initially identical, 
moving protostars. The simulated clouds that get a normal 
IMF are not uniform like this. They are always hierar- 
chically structured, or they have a hierarchy of initially 
converging velocities that lead to hierarchical fragmentation 
on the same time scale as star formation. The sink particles 
take gas from these fragments and from the condensing 
proto- fragments. 

The question we ask is whether the mass function of 
the gas reservoirs plays an important role in the final mass 
function of the stars. Even though accretion is the local 
process by which stars assemble their mass, the amount of 
accreted mass could depend more on the developing cloud 



structure (filaments, clumps, etc.) than on the presence of 
near-neighbors in a competitive environment. Everywhere 
competitive accretion occurs, cloud fragmentation also oc- 
curs at an earlier stage. Which dominates the power law part 
of the IMF? 

A test for the origin of the IMF in simulations might be 
to construct clouds with different density and velocity power 
spectra and to see if the resulting slop e of the IMF d e pend s 
on the slope of these spectra, as in the lStutzki et all (l998) 
derivation (also found in our models below). If it does, then 
cloud structure would seem to have an important role in the 
IMF. If it does not, and all density power spectra lead to the 
same IMF, then we can probably rule out cloud structural 
effects and clump mass functions for the simulated IMFs. 

One argument that cloud structure does not play 
an im portant role was given by IClark. Klessen fc Bonnell 
(2007). They suggested that if fragmentation continues 
during star formation, and if the fragmentation time is 
shorter for smaller scales, which is likely, then the IMF 
from the time-integral over all cores should be steeper 
than the instantaneous core mass function. The obser- 
vations, however, indicate tha t core mass functions have 
the same slope as the IMF (e.gjMotte, Andre fc Nerilll998l; 



Rathborne. Lada. Muench. Alves. Kainulainen fc Lombardil 
20091 ). Considering time evolution, the IMF becomes 
oc u)(k)k 2 dk, where co(k) is the formation rate of structures 
with wavenumber k, and k 2 dk is the space-fil ling partition 
function discussed above (e.g., Idi Fazio! Il986l ). For turbu- 
lence, cj(k) ~ a(k)k, where er(fc) is the turbulent speed. 
Setting cr(fc) oc k~ v for Kolmogorov turbulence exponent 
•q, we get uj(k) oc fc 1- ' 7 . If M oc fc -3 for some characteristic 
density at collapse, then we get a mass function for linear in- 
tervals of mass dn{M)/dM oc k 3 ~ v dk/dM oc M- 2 - (1 -" )/3 . 
This is the same as the mass function derived above for 
space filling structures, but now it has a slope that is steeper 
by (1 — rj) /3 ~ 2/9 for ri = 1/3 as a result of the time 
dependence. iPaloud I 20071 ) used this technique to derive 
the mass function of fragments in an expanding shell. If the 
IMF power law is not steeper by at least 2/9 compared to 
the clump mass power law, then the IM F would s e em to 
have a different origin than the clumps. iDib et al] (|2010l ) 
got around this problem by including accretion onto the 
clumps. This shifts the mass function to higher values and 
offsets the steepening from the time dependence. 



1.6 This Paper: Mass Functions in Fractal Clouds 
for Clumps that Exceed a Threshold Density 

Here we are interested more generally in mass functions for 
hierarchical and fractal clouds. We determine the mass func- 
tion for all clumps denser than a fixed value in a fractal 
cloud with a given power-law power spectrum. The depen- 
dency of the mass function on the power spectrum and den- 
sity threshold are determined. The star formation process 
itself is not considered, nor is the IMF below the plateau at 
around 0.3 Mq. Several possible connections between our 
results and the full IMF are discussed in Section [4] 

An approximately fixed d ensity threshold f or star for- 
mation makes sense physically lElmegre en (2007). Magnetic 
diffusion should suddenly become more rapid at a density 
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Figure X. (Top left)Density scan of a ID fBm cloud with two sample thresholds for clump definition. (Top right) Power spectrum of 
the density scan. (Lower left) Enlargement of a section of the fBm cloud, as indicated by the positions on the abscissa. (Lower right) 
Greater enlargement of the same cloud. 



of ~ 10 5 cm -3 because the density scaling for the elec- 
tron fraction becomes steeper, changing from n~ 1//2 to n -1 . 
This steepening occurs because charge exchange replaces 
dissociative recombination for the neutralization of ionic 
molecules, and electron recombination on neutral grains re- 
places d issociative recombination with ionic m olecules (see 
Fig. 1 in lElmegreenll979l ; lDraine fc Sutinll987l l. If the power 
of density exceeds 2/3, then the time for an initially subcrit- 
ical core t o start collapsing is comparable to th e initial free 
fall time (|Huieirat. Camenzind fc YorkelfeoOOt ). This is an 
important condition for star formation. 

Dust grains and their coupling to the magnetic field 
also change at this density. Small grains normally domi- 
nate the viscous cross section for magnetic diffusion be- 
cause of their high number density, but PAH molecules 
and small grains begin to disappear at densities in excess 
of ~ 10 cm - (|Boulanger et alJll99Ch . Gas depletion at 
high density removes ionic metals, and this also lowers the 
ionization fraction and promotes enhanc ed magnetic diffu- 
sion. The depletion time, 10 10 /n years l|Omontlll986l ). be- 
comes smaller than the dynamical time, (Gp) -1 '' 2 , when 
the density, n, exceeds 3 x 10 4 cm -3 . Grain coagulation at 



this density becomes important too, and this reduces the 
number of charged grains and the net coupling of charged 
grain s to neutrals l|Flower. Pineau Pes Forets fc Walmslevl 
l2005t ). Further coupling loss arises becaus e large grains lose 
their field line attachment at high density (|Kamava fc Nishil 
2000). All of these microscopic effects should significantly 
speed up gas collapse at a density near ~ 10 5 cm -3 be- 
cause they allow the magnetic field to leave the neutral gas 
more quickly. The mass distribution of clumps at this den- 
sity should then have an imprint on the final stellar mass 
distribution. The final mechanisms for collapse into stars 
could be less important to the IMF than the mass function 
of clumps at the threshold density where the collapse begins. 

We also study whether higher density limits produce 
steeper or narrower mass functions. Such an effect might 
explain why the power spectrum slope for the densest cores 
(which is like the Salpeter IMF slope of —2.35) is steeper 
than the power spectrum slope for low-density clouds (which 
is a ~ 1.6 — 1.8). It might also explain the apparent 
steepening of the IMF in galactic regions that have gen- 
erally low dens ities, as in low s urface brightness spirals 
l|Lee et al.l 120041 ') , dwarf galaxies ijHoversten fc G lazcbrook 
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120081; iHunter et alj|20icl). and the o uter regions of spirals 
ijMeurer et al.lll996l ; iLee et alj|2009h . The IMF would get 
steeper in these regions if stellar mass is related to clump 
mass, if there is a characteristic density at which stars form, 
and if the clump mass function gets steeper at higher relative 
densities. 

We previously found such a correlation between the 
slope of the mass function and the threshold density fo r 
clump recognition (|Elmegreenll2002l ; lElmegreen et al.| [2006). 
Both studies used three-dimensional fractal Brownian mo- 
tion (fBm) clouds with log-normal density pdfs. Here we 
experiment with fBm clouds again, first in one-dimension 
where the counting statistics are extremely good, and then 
in three dimensions, which is more relevant to star forma- 
tion. Interstellar clouds are not the same as fBm clouds, but 
the geometric resemblance is not bad if we want to address 
only simple questions. For example, fractal Brownian mo- 
tion clouds can be made with the same power spectrum as 
interstellar clouds, and they can have the s ame d ensity pdf. 
In MHD turbulence simulations, iDib et al.l (|2008h also found 
a steepening of the mass function at higher gas density, al- 
though they commented that this steepening was the result 
of stronger self-gravity at higher densities. Our experiments 
have no self-gravity and are relevant only to the geometric 
aspects of cloud structure. 

The dependence of the mass function slope a on the 
power spectrum slope P is considered in detail. Generally 
we expect steeper power spectra to correspond to shal- 
lower mass functions. This is because steeper power spec- 
tra have relatively more structure on larger scales, which 
means larger masses, so the corresponding mass function 
has relatively more high m ass clumps. Clump s are de- 
fined differently here than in IPress fc Schechten (| 19741 ) or 
Hcnncbcll e fc Chabrierl (|2008T ) . Here a clump has every part 
of it above a threshold density and the entire matter inside 
the clump boundary is counted in the clump mass. In Press 
& Schechter and in Hennebelle & Chabrier, the average den- 
sity of a clump is above the threshold but smaller regions 
inside can be below the threshold. The implications of this 
difference are discussed in section ET2"1 

iHenriksenl (| 19861 ) also suggested that the IMF comes 
from counting clumps in a fractal cloud. In his model, the 
number of clumps per logM interval is N(M) cx L~ D for 
fractal dimension D and scale L. Combined with a density- 
size relation, p cx L~ Dp , and a definition of mass, M = pL D , 
the result is N(M) cx M~ D/ ( D ~ Dp } . For constant density 
(Dp = 0), the geometric result N(M) cx M~ is recovered. 
iLarsonl (|l982l . Il99ll . 1 19921 ) also modeled the IMF for a frac- 
tal cloud but assumed M cx L for linear gas accumulations 
along filaments; then N(M) cx M~ D . Our model is similar 
to Henriksen's, but we characterize clouds by their power 
spectrum slope instead of the fractal dimension, and we con- 
sider a fixed density threshold, which makes M effectively 
proportional to L 3 . 

Larson (1973) introduced a slightly different fractal 
model in which random hierarchical cloud fragmentation 
produces a log-normal IMF. This model does not have a 
power-law component because the stellar masses are all at 
the bottom of the hierarchy, in the same level, and their 
distribution comes from a spread in the relative number of 



clumps versus interclump regions that contribute to the final 
stellar mass during the fragmentation process. Power laws 
tend to arise when only the fragments can fragment, and 
log-normals arise when both the fragments and the inter- 
fragment gas c an fragment dur ing a hierarchical process of 
star formation (lElmegreen|ll985l ). Detailed models of the first 
type were in lElmegreen! l| 19971) . The present model differs 
from both of these because we consider only the mass distri- 
bution of clumps denser than a fixed threshold, regardless of 
fragmentation history or level in the hierarchy. We assume 
no particular fragmentation model at all, only clouds with 
power-law power spectra. 



1.7 Pre-stellar Clumps 

Pre-stellar clumps have a density close to the proposed 
threshold for star formation, and they have about the same 
mass function as stars. We suggest that pre-stellar clumps 
get their mass function from fractal cloud structure as shown 
by our models here. The correspondence between pre-stellar 
clumps and individual stars is complicated and discussed in 
Section 31 There is not likely to be a one-to-one correspon- 
dence between dense clump mass and stellar mass, because 
pre-stellar clumps can fragment. Still, the full IMF, includ- 
ing the plateau and the turnover at low mass, might be ex- 
plained as a result of a combination of processes with cloud 
structure the main contributor in the power-law part. 

The mass function for lower-density clumps, like CO 
clumps in non-star-forming parts of molecular clouds, are 
not expected to have the same slope as the dense clumps 
studied here. CO requires a high column density to form 
because of molecular line shielding and dust opacity, and it 
requires a moderately high density for excitation, although 
CO line emission is possible even at lower densities if the 
CO line opacity is high. Thus CO clumps should not depend 
so uniquely on a threshold density as star-forming clumps. 
With a column density dependence in addition to a density 
dependence, CO clumps should have a shallower mass func- 
tion than clumps that are defined entirely by density. This 
is because massive low-density clumps that would not be 
included in the mass function with a pure density threshold 
can be included when there is a column density threshold in 
addition to a lower and more fuzzy density threshold. 

In what follows, Section [2] shows the size distribution 
functions for clumps in one-dimensional fBm clouds with 
various power spectrum slopes and threshold densities. Sec- 
tion [3] does the same for three-dimensional fBm clouds in 
cases with Gaussian density pdfs and lognormal density 
pdfs. A discussion of the possible relevance of these results 
and other models linking cloud structure to the IMF is in 
section 2] It is difficult to make a one-to-one connection be- 
tween clump mass functions and stellar mass functions be- 
cause of certain logical inconsistencies. An IMF that arises 
from a combination of cloud structure and late-state accre- 
tion seems to be preferred. In our model, most of the slope in 
the Salpeter IMF is from the universal geometric structure 
of turbulent clouds. 
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Figure 2. Size spectra of a ID fBm cloud for different density thresholds. The blue histograms are for the sizes of clumps exceeding the 
density threshold C and the red crosses are for the sizes of interclump regions with densities less than 1 — C. The slope is independent 
of density threshold at low size, but the maximum clump size depends on this threshold. The power spectrum slope is /3 = 5/3 for all 
thresholds. 



2 ONE-DIMENSIONAL CLOUD MODELS 

A one dimensional (ID) density strip with a pre-determined 
power spectrum was made by filling a ID strip in Fourier 
space with (0, 1) noise multiplied by k~^^ . After Fourier 
transforming this fc-space strip, the result in real space is 
a density strip with a power spectrum having a slope — f3. 
In most of the discussion below, we normalized each den- 
sity strip to have a minimum value of and a maximum 
value of 1; in some discussions, we normalize them to have 
a minimum of and an average of 0.5. To make an analogy 
with cloud mass spectra, we determined the size spectra of 
the high density regions, i.e., above a fixed threshold den- 
sity. We also determined the size spectra of the low density 
regions and of both high and low density regions together. 
The size distribution is a good representation of the mass 
distribution because most o f a cloud's mass is at a density 
close to the threshold value (ElmcErccn 2002, sec also Sect. 
[3] below). 

Figure [1] shows an example. In the top left panel is a 
random density strip 10 4 pixels long that was made with 
P — 5/3. The power spectrum of this density strip is in the 
top right panel, showing the agreement between the slope 
and the expected —5/3 slope. The bottom two panels show 
blow-up sections of the density strip. A red line indicates a 
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Figure 3. The maximum clump size versus the threshold density 
used to define a clump. The three lines are for three different total 
lengths of a ID fBm cloud. 

threshold density level of 0.5, and a green line shows the level 
0.8. We define a cloud as a region where the density is above 
a particular threshold level, so there will be one set of clouds 
with a threshold level of 0.5, and another, smaller set with a 
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Figure 4. (Top) Two sample density scans with restricted wavenumber ranges. (Lower left) The two scans enlarged and superposed 
with the same colors as in the top panels. (Lower right) The maximum clump size versus the threshold density for the blue and red 
scans, in the same colors, and for other composite wavenumbers in black and green. 



threshold of 0.8, etc.. The size of a cloud or intercloud region 
is the distance between the two points where the density 
profile meets the threshold value. Clearly there are both 
large and small clouds at all threshold levels, but there are 
far fewer clouds that stick above a high threshold level than 
above a low threshold level. 

Figure [5] shows histograms in blue of the size distribu- 
tion of clouds for various threshold levels, C, and it also 
shows histogram values, plotted as red crosses, for the size 
distribution of low-density regions for the corresponding 
threshold values 1 — C. That is, the cloud size distribution 
for C is shown in the same panel as the inter-cloud size dis- 
tribution for threshold 1 — C. Values of C are given in the 
panels. The clouds and the inter-cloud regions mirror each 
other in size distributions. This is a sensible result as the 
density strip is statistically symmetric with respect to in- 
version from top to bottom. The histograms were made for 
ID fBm density strips having 10 5 pixels, and enough strips 
were made for each C that the number of clouds or inter- 
cloud regions at a size in pixels with a logarithm between 
and 0.25 was 10 4 to maintain good counting statistics. In 
each case, /3 — 5/3. The green line in the upper mass range 



has a slope of —0.6, for comparison purposes. The rising part 
of each curve is for sub-pixel sizes. 

The maximum size of a cloud decreases with increasing 
threshold density. Higher thresholds have smaller clouds and 
larger intercloud regions. Figure [3] shows the trend of aver- 
age maximum cloud size per simulation versus the threshold 
level. The solid curve is for fBm strips of 10 5 pixels length, 
the dashed curve is for strips with 10 4 pixels, and the dot- 
ted curve is for strips with 10 3 pixels. These are for densi- 
ties that are normalized to have a minimum value of and 
a maximum of 1. For each curve, the maximum cloud size 
is essentially the whole strip length when the threshold is 
low, and the maximum cloud size decreases to zero when 
the threshold is large. The red plus symbols are for a second 
case where the density strips are normalized to a minimum 
of and an average of 0.5, with an arbitrary maximum. The 
maximum cloud sizes are about the same as in the (0, 1) 
normalization except at high thresholds, where the limitless 
maximum still produces high density regions. 

The variation in average maximum length with thresh- 
old level may be understood from Fourier analysis. Each 
density strip is the sum of sines and cosines for a wide 
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Figure 5. Mass spectra of clumps in ID fBm clouds with power spectrum slope (3 = 1 for 8 threshold levels. As in Fig. 2, the blue 
histograms are for clumps with thresholds C and the red crosses are for interclump regions with thresholds 1 — C . 



range of wavenumbers and for amplitudes that vary with 
wavenumber as k~®' 2 with f3 — 5/3 in the cases shown. 
The short wavelength oscillations modulate the long wave- 
length oscillations, and so the high intensity part of any long 
wavelength oscillation has ripples from the short wavelength 
oscillations. Figure U shows an example. The top two panels 
show the density for the full lengths of 10 4 -pixel strips with 
different Fourier components. The bottom left panel shows 
the same two strips with the same color coding magnified 
in the pixel range from 400 to 600. The bottom right panel 
plots the maximum cloud size versus the threshold level, as 
in Figure [3] for these two strips. The blue strip in the top 
and bottom left has wavenumbers only between 1000 and 
1010, i.e., it occupies a narrow range in Fourier space with 
the other Fourier components set to zero. In real space this 
strip consists of nearly pure oscillations that each go through 
the midpoint value of 0.5. The range of wavenumbers modu- 
lates the amplitude of this oscillation (top left panel) with a 
beat frequency from the difference in wavenumbers. In this 
case, the maximum width of a cloud (lower right panel, blue 
dots) is nearly independent of the threshold level for thresh- 
olds larger than 0.5 because each cloud is nearly a pure sine 
wave. Below a threshold of 0.5, the length of a cloud is the 
length of the lower envelope of beat frequency oscillation, 
which is nearly the full strip length, 10 4 . 

The red strip in the top right and bottom left of Fig- 
ure U has wavenumbers between 100 and 110 and between 



1000 and 1010. There is still the high frequency oscillation 
and modulation of this oscillation as in the blue strip, be- 
cause the same high frequencies are present, but there is 
also a low frequency modulation that pushes up and down 
the high frequency modulation (lower left panel). This low 
frequency modulation also has a varying amplitude with a 
beat frequency equal to the difference in wavenumbers (also 
10 wavenumbers). Thresholds below 0.5 (lower right panel, 
red dots) have the same long maximum length as in the blue 
strip, from the beat frequency, and above this threshold the 
short length oscillations are seen again. The maximum cloud 
length at high threshold is higher for the red strip than for 
the blue strip because the minimum frequency is lower for 
the red strip. The high threshold lengths are about equal to 
the number of pixels in the strip, 10 4 , divided by 2ty times 
the minimum k value. For the blue strip, this is 10/2-7T, or 
around 2, and for the red strip, this is 100/27T, or around 20. 

The green plus symbols in the lower right panel of Fig- 
ure U are for a case where the density strips have wavenum- 
bers only in the range of 100 — 110. The black crosses are 
for a case where the strips have wavenumbers of 100 — 110, 
310 - 320, and 1000 - 1010. These two case have about the 
same maximum cloud lengths versus threshold levels as the 
red case because the maximum cloud length at high thresh- 
old is dominated by the lowest frequencies, which are the 
same in each case. 

The lower right panel of Figure [3] resembles the lower 
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Figure 6. Mass spectra of clumps in ID fBm clouds with various power spectrum slopes, /?, and a threshold density of 0.5. The slopes 
of the histograms are indicated by the values of a — 1 and by the green lines. 



panel of Figure [3] illustrating that the maximum length de- 
creases with increasing threshold because of the isolation of 
higher frequencies at higher thresholds. The low frequency 
structure is dominated by beating effects between neighbor- 
ing frequency groups, and these effects determine the large- 
scale envelope of the density strip. This is how a cloud length 
can be greater than half the maximum wavelength in a strip. 

Cloud size distributions were also determined for other 
power spectrum indices using the density normalization 
from to 1. Figure [5] shows the case f) = 1 for a range of 
threshold densities, and Figure [6] shows histograms for many 
/?, all with a threshold density of 0.5. The slopes of the high 
mass parts of the mass functions are indicated in Figure [6] 
in the upper left corner of each panel and by the green line. 
These are slopes on the log — log plots, and therefore equal 
to ot — 1. The mass function slope gets more negative as the 
power spectrum slope gets less negative. This is sensible be- 
cause more negative mass functions correspond to relatively 



more low mass clouds, and this corresponds to relatively 
more high frequency power in Fourier space. 



Figure [7] shows the mass function slopes a versus the 
power spectrum slopes ft for the ID fBm clouds. The data fit 
the relation a = 2.83-0.66,5 between = 4/3 and 11/3. The 
threshold value was 0.5 for this. Higher thresholds would not 
change the result much because the slope a is independent 
of the threshold below the upper clu mp mass (which doe s 
dep end on the threshold). A c cordin g to lStiitzki et alj (|l998l 'l 
and iHennebelle fc Chabrierl (|2008l), the relation is a = 3 — 
/3/7 for 7 equal to the power in the mass-radius relation. In 
our ID case, 7=1 and this equation would suggest a = 
3-/3. Our actual result is close to this. 



© 0000 RAS, MNRAS 000, 000-000 



Mass Functions in Clouds and Clusters 11 




Power Spectrum Slope, (5 

Figure 7. The slope of the clump mass function versus the slope 
of the power spectrum for ID fBm clouds with a threshold for 
clump definition equal to 0.5. 



3 THREE DIMENSIONAL FBM MODELS 

Clump mass functions were also determined for three- 
dimensional clouds with and without log-normal density 
pdfs. The log-normal was made to r esemble isothermal 
clouds in simulations of MHD t urbulence I Vazq uez- Semadenil 
For 3D, we calculated clump mass in addition to vol- 
ume, although the mass turned out to be nearly equal to the 
volume multiplied by the threshold density in all cases, as 
assumed above for the ID experiments. 

The clouds were made in several steps. First we filled a 
256 3 cube in wavenumber space with random complex num- 
bers, (k x ,k y ,kz), uniformly distributed from to 1. Then 
we multiplied each number by fe _/3 ^ 2 , where k = (k 2 + k 2 . + 
k 2 ) 1 ^ 2 . We took the inverse Fourier transform of this cube 
to make a real space density distribution pm m (x,y, z). This 
distribution is a standard fBm cloud, and it has a power 
spectrum with a negative slope /3. Some of the clump mass 
functions presented below use this fBm cloud. To make a 
log-normal density pdf, we exponentiated this result: 



PLN(x,y,z) = exp(pmm(x,y,z)/po), 



(1) 



where po is for normalization. To build up the statistics of 
the cloud mass and size distributions, we generated an en- 
semble of ~ 100 fBm distributions for each of the three 
normalizations discussed below. Clumps are defined as con- 
nected pixels above a threshold density inside the overall 
cloud. The clump mass is the sum of the density values in 
all of these connected pixels, and the clump size is the cube 
root of the number of pixels. 



3.1 fBm clouds with Gaussian density pdfs 

Figure [8] shows the mass and size distributions for clumps 
in standard fBm clouds with Gaussian density pdfs in the 
case where /3 = 2.82. This f3 was chosen because it gives a 
mass function slope similar to the Salpeter IMF at a cloud- 
defining minimum density of 0.7. Three threshold densities 
are considered. The low mass and small size parts of the 
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Figure 8. Mass and size functions for 3D fBm clouds with Gaus- 
sian pdfs. The power spectrum slope is 2.82, which produces a 
mass function similar to the Salpeter IMF for low and interme- 
diate masses. The drop in the mass spectrum for high mass and 
high threshold density is a result of a superposition of low and 
high frequency wavenumbers in the fBm cloud, as discussed in 
Figure 4. Green lines are fitted to the low mass and size portions 
of the distributions. 



clump functions have about the same slope for all threshold 
densities, and there is a steepening of the functions at in- 
termediate mass and size. This steepening occurs at smaller 
masses and sizes when the threshold density is higher, as for 
the ID results. 

Figure [9] shows the slopes of the low- mass parts of the 
clump mass and size distribution functions versus the power 
spectrum slope, /?, for a density threshold of 0.7. There is 
a decreasing trend as for the ID case. Steeper power spec- 
tra produce shallower mass and size distribution functions. 
The slope of the Salpeter mass function is a = 2.35 on this 
plot, and that is associated with /3 = 2.82. This is the 3D 
in the cloud volume. The relation between the two quanti- 
ties is_O£^_4 ; 35—0 ; 715 for a threshold of 0.7. Recall that 
the lStutzki et all l |l998h and iHennebelle fc Chabrieij (|2008f ) 
result for 3D with 7 = 3isa = 3 — 0.33/3. 
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Figure 9. The slopes of the low mass and low size parts of the 
mass and size functions are shown versus the power spectrum 
slopes for 3D fBm clouds with a density threshold for clumps 
equal to 0.7. 

3.2 Comparison between our Model and the 
Press-Schechter or Hennebelle-Chabrier 
Models 

The decrease in a with increasing /3 in Figure [5] is oppo- 
site to the trend in Press & Schechter (1974) because of a 
difference in the definition of clumps. In Press & Schechter 
(1974) and in Hennebelle & Chabrier (2008), a clump is 
denned as an object with an average density above a thresh- 
old p c . The clump mass is the sum of the masses of only 
the dense parts (p > p c ) of this clump. Generally, the clump 
contains both high density (p > p c ) and low density (p < p c ) 
parts, and the high-density mass is taken equal to the to- 
tal mass multiplied by the probability that the density is 
high. This probability is the integral over the density pdf 
for all densities above the threshold. The density pdf is 
written as (2na 2 )~ 05 exp(— 0.58 2 /a 2 ) for dispersion a 2 — 
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Figure 10. (Top) Power spectrum of a 3D fBm density field made 
by the integrating the density along the z axis from z = to L 
(the size of the box in pixels). (Bottom) Power spectrum of the 
density field of a slice at Z = L/2 with one pixel thickness. The 
slope of the power spectrum of a 3D fBm density field is shown at 
the top of each panel. The projected power spectrum has about 
the same slope as the 3D power spectrum, but the slice has a 
slope that is shallower by about 1. 

(2tt)- 3 j k max P{k)dk 3 and power spectrum P(k) oc k s 

(Padmanabhan 1993). Evidently, a oc fc (3_ ' 3)/2 cc M (/3_3)/6 , 
as discussed in the Introduction. For f3 < 3, the integral for a 
depends only on the upper boundary, k max , and includes all 
k <C kmax. A.S Af increases, this upper boundary decreases 
and the range of k decreases as well. This is why the disper- 
sion of the pdf, which depends on the range of k, decreases 
with increasing M when /3 < 3. For j3 > 3, the integral for 
a depends only on the lower boundary, k min , and includes 
only k > k m in up to some cutoff in k, which may be viewed 
as the inverse of a smoothing length. As M increases, the 
lower boundary to k decreases, increasing the range of k be- 
tween that lower boundary and the upper cutoff in k. As a 
result, the dispersion a increases - opposite to the trend for 
P < 3. 

Writing either of these two k limits in terms 
of M, the clump mass function equals n(M) oc 
(M 2 cr) _1 exp(— 0.5S 2 /c 2 ), so the power law part at low M 
has a slope a = 1.5 + /3/6. This is the Press & Schechter re- 
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suit. The probability that the density is high inside a region 
increases with increasing region mass for /3 < 3 and de- 
creases with increasing region mass for /3 > 3, in proportion 
to a' 1 at low M (where the exponential is not important). 
In the first case, P < 3, this means that clumpy regions have 
dense-gas masses that progressively get smaller than in the 
purely hierarchical mass function (M~ 2 ) as M decreases. 
As a result, the clump mass function is shallower than M~ 2 
when ft < 3 in the Press & Schechter formulism. (Visualize 
an M~ 2 mass function that shifts to lower M at low M and 
higher M at high M, stretching horizontally on such a plot.) 
Similarly in the case /3 > 3, clumpy regions have progres- 
sively more mass than in the M~ 2 function as M decreases, 
making the mass function steeper (visualize the same initial 
M~ 2 function, but now compress it horizontally). This is 
opposite to the trend in Figure [9] Physically, what is hap- 
pening is that the density distribution is smoother (smaller 
<t) at larger M when ft < 3, so a region with average den- 
sity exceeding some threshold has fewer and smaller dips to 
sub-threshold density inside of it. Nearly the entire enclosed 
mass is then counted. Smaller mass regions for /3 < 3 have 
large density fluctuations inside of them, so that even when 
the average density exceeds the threshold, there are still a 
lot of regions and a lot of mass at densities less than the 
threshold. In this case, the mass fraction above threshold is 
small. 

Another difference between our model and that of Press 
& Schechter (1974) or Hennebelle & Chabrier (2008) con- 
cerns the lower limit to the mass. For /3 < 3, this lower limit 
is important, as discussed above. In the Press & Schechter 
formulism, the region defined to be a clump does not de- 
pend much on the resolution scale, which is the lower limit 
to M. Clump regions stay about the same as the resolution 
increases, and all that happens is that clump masses get 
more accurate. That is, smaller scale structure occurs near 
the density threshold and the pdf provides a more detailed 
representation of the fraction of the mass above threshold. 
In our model, the definition of each clump changes as the 
map resolution increases. More clumps and smaller clumps 
appear above the density threshold at higher resolution and 
some massive clumps disappear as they pick up fine-scale 
subregions below threshold. On the other hand, the mass 
function slope is independent of resolution in our model. 
Higher resolution extends the mass function to lower mass 
and changes the total count of clouds for fixed mass above 
the density threshold, but it leaves the mass function slope 
the same. 

Physically, our model applies when only the dense 
parts of a cloud complex are identified as clumps. These 
clumps can have highly irregular shapes, like filaments, but 
they cannot include multiple islands of threshold density 
gas separated by lower-density inter-island gas. Our model 
would identify each island as a separate clump. In Press & 
Schechter (1974) and Hennebelle & Chabrier (2008), entire 
regions with islands and inter-island gas are counted as sin- 
gle clumps if the average density of it all exceeds the thresh- 
old. In all cases, the total mass in the mass function comes 
only from the parts of the clumps that are above the density 
threshold; only the definition of a clump and the counting 
of clumps differ. Which clump definition more accurately 



applies to pre-stellar clumps depends on how observers de- 
fine their clumps and the masses of their clumps. While this 
may seem like a small detail, the slope of the clump mass 
spectrum, a, depends on the slope of the power spectrum, 
P, in opposite ways in these two cases. If clump finding al- 
gorithms change their identification of clumps as the spatial 
resolution increases, calling what used to be a single clump 
at low resolution two separate clumps at higher resolution, 
then our definition of clumps would apply. 



3.3 Projected Clouds 

Observations do not generally determine the 3D power spec- 
trum of density in a 3D cloud, but rather the 2D power 
spectrum of projected density in a 3D cloud. This difference 
between the observed and the intrinsic power spectra has to 
be investigated before we link the observed power spectrum 
to a 3D clump mass function. Two-dimensional power spec- 
tra were determined for the projected density distributions 
and for thi n slices in 3D fBm clouds. A ccording to equa- 
tion (28) in lLazarian fc P ogosvan (|2000h . the slope of the 
2D power spectrum of projected density, /3 t hick, integrated 
along a thick path through a 3D cloud, should be equal to 
the slope of the power spectrum in the original 3D cloud, p. 
The slope of the power spectrum of a thin slice in a 3D cloud, 
Ahin, should differ by 1 from the 3D slope: /3tWn — P — 1- 
The power spectrum is shallower for a thin slice than a pro- 
jection because the thin slice has more fine-scale structure. 
Figure [10] shows these 2D power spectra for a 3D fBm cloud 
with P = 2.8. The thick path projection through a 3D fBm 
cloud (top panel) was determined by integrating along the 
z axis from z = 1 to z = L, which is the 3D box size. The 
thin slice projection (bottom panel) was determined by us- 
ing only the z — L/2 position to make the 2D map. Figure fTUl 
indicates that the 2D power spectra have approximately the 
slopes that are expected. The projected density field has a 
2D power spectrum slope /3 t hick = —2.76, and the 2D power 
spectrum of the thin slice has a slope /3thin = 1.92. These 
are approximately equal to, and 1 less than, respectively, the 
3D slope of P = 2.8. 

This result for 2D indicates that observation of 2D 
power spectra for projected or line-of-sight integrated clouds 
having slopes of about —2.8 (Table 1) correspond to 3D 
power spectra inside those clouds with about the same slope. 
Correspondingly, Figures [8] and [9] now suggest that for this 
slope, the mass function of clumps has a slope equal to the 
Salpeter value, a = 2.35. Thus the power-law portion of the 
stellar IMF can have its fundamental origin in the geometric 
structure of clouds. Section [4] discusses the implications of 
this point in more detail. 

Mass functions for clumps observed in projection are 
shown in Figure [TT1 These were made from the same 3D fBm 
clouds used for Figure [9] but the search algorithm was done 
on 2D projections of the clouds made by integrating along 
the z axis. Mass functions for the 3D fBm clumps shown 
previously are reproduced as dashed lines. Three different 
density thresholds are considered with a 3D power spectrum 
slope p — 11/3. This is the slope of the 3D power spectrum 
for velocity in classical 3D incompressible Kolmogorov tur- 
bulence. The green lines in Figure [TTJ show fits for a density 
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Figure 11. Projected mass functions of clumps in 3D fBm clouds 
obtained by integrating along the z axis. The mass spectra are 
shown by solid lines for thresholds of 0.7, 0.8 and 0.9. The mass 
functions of the corresponding 3D fBm clouds are shown by 
dashed lines. The power spectrum slope is /3 = 11/3. Green lines 
are fitted to the 2D and 3D mass functions for a density threshold 
of 0.7. 

threshold of 0.7. As for the 3D mass functions, the slopes of 
the low-mass parts of the 2D mass functions are about the 
same for all density thresholds; the slopes of the upper mass 
parts are steeper, and the upper mass limits are smaller, for 
higher density thresholds. The slope of the 3D mass func- 
tion is written in the Figure as a = 1.9. This is also the 
value expected from Figure|5]for /3 = 11/3. The slope of the 
projected 2D mass function is «2D = 1-4. Evidently, the 2D 
mass function of the projected 3D cloud is shallower than 
the 3D mass function in the unprojected cloud by ~ 0.5 in 
the slope. This implies that the 2D mass function has pro- 
portionally more high mass clumps, presumably from the 
blending of smaller clumps. 

The 2D mass function of projected clumps is not nec- 
essarily the same as the mass function of interstellar clouds 
observed with dust emission integrated along a long line of 
sight. If there is blending on the line of sight because each 
clump is relatively large in angle, as would result from a low 
threshold density for the definition of a clump, then this pro- 
jection result might apply. However, if the density threshold 
is very high and there is no blending for a cloud with a finite 
line-of-sight path length, then the observed projected mass 
function could be closer to the true 3D mass function. 



3.4 fBm clouds with log-normal density pdf's 

A second set of models was made with log-normal density 
pdfs, which are exponentials of the previously discussed fBm 
clouds (Equation [1] see also lElmegreenl (|2002T l). We consid- 
ered three types of density normalization. The first forced 
each log-normal cloud to have a minimum density, pln, of 
and a maximum density of e 1 . This was done by selecting po 
in equation {TJ to be the maximum value of PfBm, designated 
as pfBm.max- Each fBm cloud had a different po = pfBm,max- 



Clumps were defined to be regions where the density pln 
was larger than some fixed fraction, /dump, of the peak value 
of e 1 . The same fraction was used for all of the log-normal 
distributions in the ensemble; different / c i ump were used to 
see how the mass distribution depended on the threshold 
density. 

A second case of pln normalization had a minimum 
density of and a different maximum density for each cloud. 
The normalization po was now taken to be the same for 
each cloud and equal to the average of all the pfBm.max 
from the 100 previous ensembles, po =< ptBm,max >■ Thus 
the density maximum in a log-normal cloud was pln, max = 
exp(p fBm ,max/ < PfBm,max >)■ The random numbers were 
the same as in the first case. The density threshold for 
the definition of a clump was set equal to a fraction /dump 
of each density peak, i.e., / c iumpPLN,max. Because the den- 
sity peaks are all different in this normalization, the density 
thresholds are all different too, although the density thresh- 
old fractions are the same. The fractions are taken to be 
the same as in the first normalization case, with a different 
/dump for each ensemble of 100 clouds. 

The third case of pln normalization had a minimum 
density of and a variable maximum pln. max defined in the 
same way as in the second case, i.e., with the same po =< 
Pffl m ,max > for each cloud. Now the density threshold to 
define a clump was taken to be the same for each cloud, not 
variable with the same fraction of the peak. 

The motivation for considering these different cases is 
our hypothesis that cloud collapse to stars begins to get 
significant at a certain ch aracteristic densit y where the mi- 
croscale processes change (|Elmegreenll2007l ). We would like 
to know, for example, how the clump mass spectrum de- 
pends on the peak gas density in a kpc-size region for a fixed 
characteristic density (this is closest to case 3). We would 
also like to know how it depends on the relative density of 
the collapse threshold (close to case 2). If the absolute den- 
sity for the onset of collapse is the same everywhere, then 
the timescale for collapse would be about the same and the 
star formation rate should scale as the first power of the av- 
erage ISM density. If, on the other hand, the relative density 
for collapse is the same everywhere, then the timescale will 
scale with the inverse square root of the average density and 
the star formation rate will scale with the 1.5 power of the 
average density. 

Figure [12] shows the mass and size distributions and the 
mass-size relation for the three different cases using f3 = 5/3. 
Each curve in the top two rows is a different clipping level, 
with lower levels producing larger numbers of clumps. The 
clipping levels were chosen to avoid excessive cloud merging 
at low levels and too few clouds at high levels. Each dot 
in the lower row represents a single clump. Lower clipping 
levels produce lower mass clouds for the same size. 

The mass and size distributions are related to each 
other by the mass-size correlation. For mass distribu- 
tion dn(M)/dM = M~ a , size distribution d\(R)/dR = 
i? _A , and mass-size relation M oc R"' , we should have 
dn{M)/dM = (dA(R)/dR)(dR/dM), from which we get 
a = 1+ (A — l)/7. The bottom panels of Figure [12] indi- 
cate that 7 ~ 3.0 for all cases, which implies that the average 
density of clouds is nearly uniform when there is a threshold 
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Figure 12. (Top) Clump mass spectra of 3D fBm clouds with log-normal density pdfs and power law slopes = 2.8. Three threshold 
values equal to 0.6, 0.7 and 0.8 are considered. (Middle) Size distributions of the clumps. (Bottom) The mass-size relations. The three 
methods for cloud normalization, corresponding to the three columns, are described briefly in each panel and discussed in the text. 



density to define the clouds. This average density is approx- 
imately 1.08 times the threshold density, depending a little 
on p, because most the cloud volume is near the cloud edge 
where the density equals this threshold. Recall that clumps 
made from fBm clouds are fractals with convoluted surfaces 
near the threshold density, so a higher fraction of the mass 
is near this threshold density than in an isothermal sphere. 
For 7 = 3.0, a = 1 + 0.33(A - 1). All of our 3D results 
approximately follow this expression. 

The mass functions and size functions in 3D resemble 
those in ID and in the regular fBm case in the sense that 
they have similar slopes at low mass but then steepen to dif- 
ferent degrees at high mass. The mass where the steepening 
begins is lower for higher threshold densities. This trend is 



present for all three cases in Figure [121 The average slope of 
the mass function is a blend of the common low-density slope 
and the steepened high-density slope. This average steepens 
from —2.1 to —2.8 as the relative threshold density increases 
from 0.6 to 0.8 in the first normalization case. The average 
slope steepens in the second case from —2.12 to —2.76 and 
in the third case from —1.94 to —2.27. 

The steepening with increasing threshold density is low- 
est in the third case, which is the one with a variable peak 
and a constant threshold for clump definition. The variable 
relative density threshold that results in this case effectively 
blends together the mass functions from cases with higher 
and lower absolute thresholds. 
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4 DISCUSSION: HOW IMPORTANT IS 
CLOUD STRUCTURE TO THE IMF? 

As a step toward understanding the IMF, we have investi- 
gated the mass functions of clumps that exceed a threshold 
density in fractal Brownian motion clouds. We expected that 
mass functions will resemble volume functions when there is 
a threshold density, and we verified that here. We also found 
relations in ID and 3D between the slopes of the power laws 
for the mass functions and the power spectra, reminiscent of 
the analytical results in Stutzki, et al. (1998) and Hennebelle 
& Chabrier (2008), although not exactly the same. For the 
most realistic case, a 3D fBm cloud with a projected power 
spectrum similar to what is observed in molecular clouds, 
the mass function of clumps denser than a threshold is a 
power law with the Salpeter slope. This result does not de- 
pend on the threshold density much, nor on the projection. 
Only the maximum clump mass depends on the threshold 
density. As a result of this agreement between the clump 
mass spectrum and the stellar IMF for realistic structure in 
interstellar clouds, we are led to believe that the power-law 
part of the IMF results in large part from cloud structure 
at a threshold density. The low-mass part of the IMF would 
then arise from an inefficiency of star formation below a cer- 
tain mass, like the thermal Jeans mass. We consider this in 
more detail here. 

The clump function is a reasonable starting point for 
studies of the IMF, but the IMF cannot be explained 
so simply. If there is a fixed dens ity threshold (e.g., 
Ijohnstone et all 1 12004 lElmegreedl2007h , then larger clumps 
with the same temperature contain more Jeans masses and 
should fragment into more sub clumps and stars. The as- 
sumed one-to-one correspondence between clump mass and 
star mass is lost. This is a problem faced by observations 
of the pipe nebula |kathborrie ct al. 200^): the clump mass 
distribution resembles the IMF, but the largest clump con- 
tains many stars and breaks the assumption of a one-to-one 
correspondence between clump mass and star mass. On the 
other hand, if there is such a correspondence, i.e., if the Jeans 
mass scales with the clump mass, then massive stars form 
in lower density clumps or in the low-density periphery of 
clumps (in the absence of temperature variations) , and this 
positioning is contrary to the observed mass segregation in 
young clusters. More physical processes are required, such 
as heating. 

An additional problem is that if larger clumps contain 
more Jeans masses and form more stars, then there would 
be a stellar mass function inside each clump, which could, 
in principle, differ from the clump mass function. Only the 
sum of the IMFs inside all of the clumps has to equal the 
total cloud IMF, which is the generally observed function. 
There is virtually no systematic study of individual (local) 
IMFs inside embedded-star clumps. 

There are several ways in which the sum of individual 
clump IMFs can give the observed IMF. One is for the IMF 
inside each clump (the "internal clump IMF") to be the 
same as the cluster IMF, up to a certain maximum star 
mass dependent on the clump mass, and for the clump mass 
function to have a slope shallower than or equal to a = 2 
llElmegreen. Elmegreen. Chandar. Whitmore fc Reganl 
2006). This seems ruled out by the observation of clump 



mass functions with steeper slopes, similar to that of the 
IMF. A steeper clump function produces a summed IMF 
that is steeper than the individual clump IMF (see below). 

Another way to give the right summed IMF is for the 
internal clump IMF to be a random sample of the cluster 
IMF up to arbitrarily high stellar mass in each clump. Then 
the clump mass function does not enter into the IMF. How- 
ever, it is unrealistic to expect each clump to be able to 
produce a star of arbitrarily large mass. With additional 
physical processes, such as continued accretion from other 
clumps, random IMF sampling might be possible, but not 
with static clumps as the only reservoirs for internal clump 
IMFs. 

In a third way, the internal clump IMF can be an ar- 
bitrary function of the star-to-clump mass ratio. Then the 
summed IMF will have the same power as the clump mass 
function in the power-law range, regardless of the internal 
clump IMF. In this third case, the universality of the IMF 
follows from the universality of cloud structure. The present 
study of fBm clump mass functions fits closest to this third 
scenario. There need not be a one-to-one correspondence be- 
tween clump mass and star mass. Instead, there has to be an 
internal IMF for each clump that is relative, i.e., the stellar 
masses scale up or down with the clump mass. This might 
be the case if competitive accretion dominates the assembly 
of stars inside each clump, and if heating, merging, ejection, 
and other processes contribute to the final stellar IMF inside 
each clump, regardless of what that IMF is, as long as all of 
the rates and masses are proportional to the clump mass. 

These three scenarios are worth reviewing here. We con- 
sider in all cases a constant ratio of total star mass to clump 
mass in each clump, so the total star mass in each clump 
is equal to some fraction e of the clump mass. We also con- 
sider for simplicity only the power-law part of the IMF, ig- 
noring the low-mass plateau. The clump mass function is 
dn c (M c )/dM c = n c oM~ a , as before. 

For the first of the three scenarios, we consider a 
star mass function inside each clump of dn 3 (M s ) /dM s = 
n s oM^' 1 ~ x for Salpeter x = 1.35. The integral of M s times 
the star mass function from a fixed minimum stellar mass 
Afs.min up to some very large stellar mass has to equal the 
total stellar mass in the clump, which is eM c . That means 
n s0 = AeM c for A = (x — i)M*^ in , and it means the IMF 
in a clump of mass M c is n a (M a \M c ) = AeM c M s ~ 1 ~ x . The 
summed star mass function for all of the clumps is the inte- 
gral of M s n s (M S \M C ) over the clump mass function from the 
minimum clump mass that is likely to contain a star of mass 
M s , which is M c , m in = xM£/(Ae), up to the maximum 
clump mass in the cloud, M c , max . The minimum clump mass 
comes from the requirement that J M n B oM~ 1 ~ x dM s = 1. 
The summed star mass is therefore 
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for For a = 2, it is 

n S u m (M s ) = AericoM- 1 -* In (AeM c>max / [xMf]) . (3) 
When a sj 2, the first term in the parenthesis of equa- 
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tion ([2} dominates the second and the summed stellar mass 
function is oc M~ ~ x , which is the same as the individ- 
ual clump IMF. If a = 2, the same is true approximately, 
with only a logarithmic deviation. If a > 2, then the second 
term in the parentheses of equation ((2} dominates the first, 
and the summed IMF is oc M~ 1 ~ x(a ~ 1 \ This summed IMF 
can be much steeper than the Salpeter function if a > 2 
jKroupa fc Weidned I2003T ) . Because dense clumps are ob- 
served to have a > 2, individual clump IMFs have to be 
shallower than the whole-cluster IMF, which is the Salpeter 
function. 

For the second of the three scenarios, a star of any mass 
is supposed to be able to form in a clump of any mass, in 
which case M C:m in is constant in the integral of equation 
[21 independent of M. Then the summed IMF is trivially 
proportional to the individual clump IMF, M s ~ 1 ~ x . 

The third scenario allows a variety of stellar mass func- 
tions inside each clump and always gives a summed stellar 
function with a slope that is equal to the slope of the clump 
mass function. To demonstrate this, we assume the internal 
clump IMF is a distribution function, P(f)df, of the ratio 
of star mass to clump mass, / = M s /M c for / between f m in 
and f m ax- For simplicity, P(f) = Pof~ p - Then the summed 
stellar mass function is 

n SU m(M B ) = J P(f)n c (M c )df 

P Ml~ p Mr a ~ 2 dM c , (4) 



which integrates to 

.... P M}~ P 
n S um(M s ) = 



— (M p ~ a ~ 1 

j\ \ c,ma 



(p - a - 1) 



M 



p—a— 1 \ 



(•>) 



If clumps larger than a certain mass Mj fragment into 
stars, then M c , min = MAX(Mj, M 3 /f max ). Also, M c , max = 
M s f m 

in • 

Now consider clump masses larger than the fragmenta- 
tion limit, i.e. M 3 /f max > Mj. Then M c , min = M a /f max , 
and the stellar mass function is 
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This IMF is proportional to M^ a , which is the same as the 
clump mass function with a slope that is independent of 
the individual clump IMF slope p. Below the fragmentation 
limit but not as low in stellar mass as the minimum fragment 
size, we have M Cjmin = Mj and also M s > Mjf min . Then 
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where the second term in the brackets is less than the first 
for 1 + a > p. In this case, the summed IMF is proportional 
to M]~ v , which is the function P(p) in intervals of log P. 

This result for the third scenario implies that if clumps 
need to contain at least one Jeans mass before they fragment 
into a little sub-cluster, and if the Jeans mass, Mj, is about 
the same for all clumps (i.e., for a constant temperature and 



threshold density), then the summed IMF above Mjf ma x, 
which is the largest star that comes from an Mj piece, is a 
power law with the same slope as the clump mass function. 
Below Mjfmax, the summed IMF is the internal clump mass 
function. If p ~ 2, then the IMF plateau may be explained 
as the part below the smallest unstable clump. It would 
be composed of stars coming from that smallest unstable 
clump, and from other clumps with slightly larger masses 
(|Elmegreenll200d ). 

What type of model will produce a relative stellar mass 
function, p(f), inside each clump, rather than an absolute 
stellar mass function inside each clump? Such a function 
would seem to result from self-similar processes that scale 
to the clump mass, and that would include protostellar ac- 
cretion of a gas mass proportional to the clump mass, and 
protostellar interactions, including possibly coalescence. As 
all of this occurs in clumps more massive than Mj, the mo- 
tions would be transonic and Mj itself would be uncon- 
nected with the final star mass. Note that in this model, 
the number of fragments in a clump would not be propor- 
tional to the number of Mj in the clump, because that would 
give stellar masses that are centered around Mj rather than 
stellar masses that scale with the clump mass. In this inter- 
pretation, Mj enters the final IMF in its establishment of 
a minimum fragmentable clump mass at the IMF plateau. 
Also in this model, massive clumps would not produce more 
stars than low-mass clumps with the same average stellar 
mass, but they would produce the same number of stars as 
the low-mass clumps and a higher mass for each star. Thus 
the third model is like assuming each clump produces stars 
with masses proportional to the clump mass. It is effectively 
the same as assuming a one-to-one correspondence between 
clump mass and star mass, even though each clump forms 
several stars with a local, and even arbitrary, relative IMF. 



5 CONCLUSIONS 

Interstellar clouds have a power spectrum of column den- 
sity that is a power law with a negative slope in the range 
from ~ 2.7 to 3. The underlying structure is not a clump- 
interclump medium but a continuum of densities with struc- 
tures in the form of sheets, filaments, and clumps. Surveys 
that are sensitive to column density or density view these 
clouds as collections of discrete clumps, which are the dens- 
est parts of the continuum. In many models, star formation 
also selects the densest parts of a cloud. Here we investigated 
a model for the IMF based on this principle. We determined 
clump mass spectra above threshold densities for clouds with 
power-law power spectra. The results showed a dependence 
of the clump mass spectrum on the slope of the power spec- 
trum, as expected from analytic theory, and a trend toward 
decreasing maximum mass with increasing relative threshold 
density. For realistic power spectra, an IMF-like slope results 
for the derived clump spectra. At a very basic level, this re- 
sult supports numerous models of star formation where the 
intermediate-to-high mass part of the IMF has its basic form 
from cloud geometry. 

The IMF has to be more complicated than this, how- 
ever. If there is a characteristic density for star formation 
and a wide range of clump masses above that density, as is 



© 0000 RAS, MNRAS 000, 000-000 



18 M. Shadmehri and B. G. Elmegreen 



the case for clouds with power-law power spectra, then larger 
clumps contain more unstable sub-clumps and should form 
more stars or more massive stars. If they just form more 
stars with the same characteristic mass, e.g., the thermal 
Jeans mass, then the assumed connection between clump 
mass and star mass is lost, as is the connection between the 
clump mass function and the IMF. The IMF needs a differ- 
ent origin in that case. If, on the other hand, the average 
star mass inside a clump scales with the clump mass, then 
the clump mass spectrum can be the same as the IMF in 
the power-law part. Such scaling allows for the interesting 
possibility that clumps have internal IMFs that are differ- 
ent from observed IMFs in whole clusters, which is the sum 
from many clumps. Internal IMFs inside clumps have not 
been investigated much so this possibility remains open. 

Competitive accretion and clump coalescence play an 
important role in forming the IMF, as shown by many nu- 
merical simulations. These processes usually operate along 
with cloud fragmentation, so it is unclear whether the bulk 
of the IMF power law still comes from geometry. That is, 
the basic dn/dM oc M~ 2 power law that arises in so many 
fragmentation models is already very close to the IMF (the 
Salpeter slope is = —2.35 in this convention). Accretion, 
stellar winds, coalescence, and other microphysics need only 
modify this initial function by a small amount, giving a slight 
bias toward lower mass objects. Experiments of competitive 
accretion in clouds with different density power spectra or in 
uniform clouds might reveal the overall importance of cloud 
geometry in the IMF. 
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